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We study the quenched and unquenched lattice Schwinger model with Wilson fermions. The lowest non-trivial 
order of the systematic expansion recently proposed by Sexton and Weingarten is shown to allow good estimates 
of long distance physics from quenched configurations. Results for the static potential and the lowest bound state 
mass are presented. 



1. INTRODUCTION 

Renewed efforts have recently been made to de- 
sign improved algorithms for the simulation of 
theories with dynamical fermions. The pressing 
nature of the problem in QCD is well known. 
Since much of the effect of dynamical quarks in 
QCD is predicted by perturbative arguments to 
amount to a shift in the effective lattice spacing, 
it will require considerable computational power 
and improved algorithms to isolate meaningful 
and non-trivial effects for quantities of physi- 
cal interest. The Schwinger model offers a suit- 
able testing ground of these questions for sev- 
eral reasons. The quenched lattice theory has 
(unlike QCD) an almost trivial pure gauge sec- 
tor with most of the dynamics contributed by 
the fermions. The effects of unquenching are ex- 
pected to be pronounced. Much is known exactly 
and perturbatively about the corresponding con- 
tinuum theory. Obviously the Schwinger lattice 
model is computationally more accessible than 
full QCD. 

Sexton and Weingarten have recently proposed 
a systematic expansion of the fermion determi- 
nant in terms of its constituent gauge invariant 
loops [1],[2]. An initial test with QCD on rel- 
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atively small systems (6 4 ) has given promising 
results in comparison with hybrid Monte Carlo 
simulation of the full theory. They showed [2] 
that the discrepancy between Wilson loops mea- 
sured in quenched and full simulations can be 
accounted for by the lowest non-trivial system- 
atic correction estimate given by this expansion. 
This might seem surprising at first sight. The 
systematic correction, or deficit, is estimated at 
lowest order from an approximate evaluation of 
the overlaps of an observable with the trace log 
of the fermion matrix and with the smallest Wil- 
son loop (plaquette) - both with respect to the 
quenched vacuum. One might expect this esti- 
mate to be severely tested when truly long range 
physics is studied. With this in mind, we have be- 
gun a study of the Schwinger model. We wish to 
test the limits of the method and to uncover any 
physical insights relevant to the Schwinger model 
itself and to the continuing work on QCD. 

2. METHOD 

Details of the method can be found in [1],[2]. 
The implementation for the Schwinger model is 
essentially the same as for QCD. We briefly re- 
view some of the basic ideas and formulae. We use 
the Wilson formulation of the lattice Schwinger 
action 

S = S G + S f where S f = ~^M[U]ip (1) 
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in a suppressed index notation. The fermion 
matrix M can be decomposed into its red/black 
(even/odd) components 

1 kK° 
nK eo 1 



M = 1 + kK 



(2) 



From this one may construct a Hermitian, posi- 
tive definite, red/black preconditioned matrix 

H (e) = M (e) M (e)t j M (e) = \ - K ? K e ° K° e (3) 

from which one may reconstruct the effective 
fermion-gauge action for nj flavours of fermion 
5 e fj = —rifT/2, where 

T = Trln(MMt) = 2 In det M = In det H e . (4) 

We may formally expand T in terms of a Gram- 
Schmidt orthogonalised basis of ordered gauge- 
invariant loops, suitably symmetrised over spatial 
orientations and translations: 



T = ^2 = a t S t 



R„ 



(5) 



t=o 



where L n represents a truncation after n terms 
and R n the corresponding remainder. Here So = 
1, (loops of size 0) and 

S 1 =S 1 -<S 1 > (6) 

where Si = Sg is the (unnormalised) plaquette 
action. The higher loops S n used to form the 
orthogonal basis are also unnormalised sums (e.g. 
S2 on a unit L x L lattice has value 2L 2 ). The 
basis vectors themselves satisfy 

< SiSj >=< Sf > 6ij , < Si >= (i > 1) . (7) 

It is immediately obvious that the truncation Lo 
is the quenched approximation and that L\ cor- 
responds to a shift in gauge coupling 6/3. Each 
truncation gives rise to an effective partition func- 
tion 



Z n = l[dU]exp(-S G + ^-L n ) 



(8) 



and the corresponding expectation value < A > n 
of an arbitrary operator A. One can show [1],[2] 
that the deficit in such an estimate due to the 
neglect of R n is approximately 

A < A > n = Ul ' 



-(< AR n > n 
< A > n < R n >„) 



In view of the above arguments, < A >i can be 
estimated from quenched configurations shifted 
in [3. The corresponding deficit can then be ob- 
tained from 



A < A >i 



(< AT >i - < A >i< T >i 



ai < ASi >i) 



where 



ai 



26£ _ < T5i >i 
n f ~ < S? >i . 



(10) 



(11) 



(9) 



The stochastic evaluation of T = TrlnM using 
Gaussian noise vectors and a Chebychev approxi- 
mation to the logarithm is described in [1]. In the 
present work, we find 40 noise vectors are neces- 
sary to obtain fluctuations which are comparable 
with or smaller than those with respect to the 
gauge fluctuations. 

An immediate test of the method is to compare 
< A > from a full dynamical fermion simulation 
(FDFS) at /3 f = 13 +613 with < A > 1 +A < A > 1 
estimated from quenched configurations at /?o [1]. 
For this purpose, we take A in turn to be Wilson 
loops, the static inter-fermion potential V(r) de- 
rived from these, and m v the mass of the 'vector' 
bound state. The latter is known in the massless 
continuum model to be given by ef^/v where e is 
the fermion charge. On the lattice we deduce it 
from the large Euclidean time dependence of zero 
momentum correlators involving ipli'4'- These 
were constructed from 'quark' propagators which 
were obtained by a conjugate gradient-style solver 
applied to the red/black preconditioned matrix 
if( e ). We tried further preconditioning in mo- 
mentum space with a suitably tuned free fermion 
propagator. For the present application the ac- 
celeration achieved did not repay the overheads. 

The FDFS was performed with the hybrid 
Monte Carlo algorithm [3] which is exact for even 
nf. We choose nj =2 for the comparisons. The 
continuum version of the massive multi-flavour 
model is known to have a non-trivial spectrum 
on finite systems with light fermions [4]. In par- 
ticular the limits L — ► 00 and m —> do not com- 
mute. This will be a challenge for future detailed 
analyses of the flavour dependence of unquenched 
models. 
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3. RESULTS 

Our preliminary analysis was performed on lat- 
tices of size 32 2 , at (3 values from 2.0 to 3.0, and 
k from 0.25 to 0.265. For example, we generated 
2000 quenched configurations at [3 = 2.5 which 
we first gauge-fixed (Landau) in case Fourier ac- 
celeration was required later. We used a Lanczos 
algorithm to estimate the maximum and mini- 
mum eigenvalues of at each k (needed for 
the Chebychev approximation) and then accumu- 
lated stochastic estimates of T and T 2 . We also 
accumulated measurements of the operators A on 
each configuration. From these, we evaluated 6/3 
for 2 flavours (Eqn. 11), some higher expansion 
coefficients a n (such as two and three plaquette 
operators) and deficit estimates A < A >i for 
the above observables. 

At /3 = 2.50 and k = 0.26, we found 6/3 = 
—0.212 ± 0.010. These and other all other errors 
were evaluated using bootstrap with subensem- 
bling to detect and remove autocorrelations. We 
then generated a further 2000 configurations in- 
corporating dynamical fermions (k = 0.26) at 
f3 = /?o ± 8/3 = 2.29. A simple preliminary 
test is to compare < Si > (=< Sg >) for the 
quenched data, 783.3 ± .3 with the full fermion 
data, 782.8 ± .3 at shifted [3. Since the above 
uncertainty in 6/3 translates to one of ±1.4 in 
< Sq > this is a very satisfactory comparison. 

Fig. 1 shows the vector state effective mass 
m v obtained from the FDFS at [3 = 2.29, k = 
0.26 in comparison with that estimated from the 
quenched data together with the deficit (10). The 
dashed line in this figure serves both to indi- 
cate the quenched value at [3 = 2.29, and the 
Euclidean time range used to extract mass esti- 
mates. From the FDFS and from the corrected 
estimate from the quenched data respectively we 
find m v = 0.369(3) and 0.367(7) where the brack- 
eted digit is the error in the last significant digit. 
These are in complete agreement. In the cor- 
rected value, the contribution due to the [3 shift 
is —.041 and that due to the truncation deficit, 
-.037. 

Fig. 2 shows a comparison of the static po- 
tential using quenched data at [3 = 2.50, FDFS 
data at 2.29 and the corrected estimate from the 
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Figure 1. Effective masses from the FDFS and 
corrected quenched data at k = 0.26. 



quenched data. The results have been expressed 
in dimensionless units using our result m v = 0.37 
in lattice units. Also shown is the (exact) in- 
finite volume quenched result (dashed line) and 
the continuum single flavour massless Schwinger 
result (full line) obtained by evaluating Wilson 
loops due to a vector field of mass e/^/ir. The 
continuum result is 
V{m v R) 



,-(1 



-m v R 



(12) 



The comparison is satisfactory although the er- 
rors associated with the deficit estimate (10) get 
quite large at large distances. There could be a 
systematic under-estimate at increasing R, but 
this is not yet statistically significant. The same 
effect is seen in a direct comparison using Wil- 
son loops with increasing area (not shown here). 
However, it is clear that the potential screening 
effects due to light dynamical fermions are capa- 
ble of being reproduced by these low order esti- 
mates. 

In Table 1, we collect information on the rel- 
ative contributions to the expansion of T = 
2TrlnM. We display the coefficients 01,02, as a 
and 03j (3a corresponds to three plaquettes in 
a row, 3& is a 'chair' of three plaquettes). Also 
listed is the corresponding root variance a n of S n 
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Figure 2. Static potential in physical units: 
quenched at [3 = 2.50 (crosses), corrected for the 
deficit at first order (open circles) and from FDFS 
(full circles). The curves are described in the text. 



ably large. Having estimated the (small) higher 
coefficients in the expansion, one may now con- 
sider simulating with a higher order effective ac- 
tion. This could be done in the same spirit as 
above i.e. by making use of a deficit correction 
or, perhaps more promisingly, by using the effec- 
tive action to generate candidate configurations 
within an exact simulation. These and other is- 
sues will be addressed in future work. 
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Table 1 



Contributions to the fluctuations in T = 2TrlnM. 







Cn 


Product 


Si 


0.212(10) 


10.7(2) 


2.27(11) 


s 2 


0.085(7) 


15.2(2) 


1.29(11) 


S3a 


0.041(6) 


14.2(3) 


0.58(9) 


S3b 


0.025(4) 


25.2(4) 


0.63(10) 



and the product a n a n . It is clear that the sin- 
gle plaquette provides the largest contribution to 
the fluctuations in T and that the coefficients are 
decreasing. It is presumably for this reason that 
we can obtain reliable estimates of the remain- 
ing deficit incurred when working with quenched 
configurations at shifted (3. 

4. COMMENTS AND CONCLUSIONS 

The Schwinger model results reinforce those al- 
ready obtained in QCD [1,2]. One may indeed 
use quenched configurations to study non-trivial 
physics due to dynamical fermions. However, 
the work required in the stochastic evaluation of 
TrlnM using current techniques is still unaccept- 



